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Interaction effects in multi-subband quantum wires 
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We investigate the effect of electron-electron interactions on the transport properties of dis- 
ordered quasi one-dimensional quantum wires with two or more subbands occupied. We apply 
two alternative methods to solve the logarithmic divergent problem, namely the parquet graph 
theory and a renormalization group calculation. We solve the group equations analytically in 
the weak coupling limit and find a power-law for the temperature dependent conductivity of 
a multi-channel system. The exponent is roughly equal to the inverse of the number of the 
occupied subbands. 
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One-dimensional interacting electron systems have at- 
tracted considerable interest recently. Much attention 
has focussed on striking effects in presence of disorder 
potentials such as a vanishing conductance at zero tem- 
perature and non-Luttinger liquid behaviorEP. Experi- 
mentally, the temperature dependence of the conductiv- 
ity <j(T) of long wirescf seems to favor predictions of the 
Tomonaga-Luttinger mode©. Although the latter has 
become a standard way to describe correlated Id elec- 
trons, there are still arguments for ordinary Fermi liq- 
uid theory to be valid in the presence of disordeitf. In- 
vestigations of multi-component models where directed 
mainly towards a small number of coupled chainscP , trig- 
gered by theupossible relation to 2d superconductivity® , 
or tunnelingEP . On the other hand, in transport exper- 
iments one is able to drive a system from the extreme 
quantum limit of one to a large number of TV occupied 
subbands. The question then is how Coulomb effects at 
low temperatures change when TV is varied, with the pos- 
sible re- establishing of Fermi- liquid behavior in the limit 
of large N. 

In this Letter, we present a microscopic theory for the 
conductivity a of a multi- subband wire in presence of 
an arbitrary, weak static disorder potential and Coulomb 
interactions between the electrons. One of our central 
predictions is a temperature dependence 



a (T) ~ T^ N 



(1) 



in the limit of weak interactions a*/W 3> 1, where a 
is the effective Bohr radius and W the wire width, 
similar 1/N dependence has been found by Flensbergl 
for the conductance of a quantum dot coupled to quan- 
tum point contacts in the framework of a multichannel 
Tomonaga-Luttinger model. The result Eq. ([]]) is a gen- 
eralization of the one-cannel casetf cr(T) ~ T x ( 2g2 ~ 9l \ 
where the factor A<?2 ~ 1/2 turned out to be 'nearly' 
universal, while Xgi — c\/a*kp (c\ a constant of order 
unity which depends on the structure of the wire) is neg- 
ligeable since typically kp ~ l/W. n n 

We generalize the parquet diagram theoryErtP to the 



present case and give an alternative, microscopic deriva- 
tion using a standard one- loop renormalization group 
(RG) calculation which reproduces exactly the vertex 
equations of the parquet theory. It turns out that the 
singular (q — > 0), logarithmically divergent Coulomb in- 
teraction has to be renormalized. Namely, one has to 
take into account all diagrams which contain powers of 
the interaction matrix elements, which physically means 
that the interaction effectively is replaced by a screened 
Coulomb interaction. This reduces the problem to that 
of short range interactions which dominate the low-T 
behavior of the conductivity. The latter conclusion is 
in contrast to calculationsEfP , where a contribution from 
4fcp scattering was found to be dominant at low T in 
the case of long range Coulomb interactions and delta 
impurities. At present, it is not clear if this 'Wigner- 
crystal pinning' effect dominates over the 2fcF-pContribu- 
tion in experiments with relatively short wires© . Never- 
theless, we believe that there should be a wide tempera- 
ture regime over which the 1 /N behavior Eq. (|l]) can be 
observed. 

We start from a multi-band quantum wire with dis- 
persion E n k, where n is the band index and k the mo- 
mentum of plane waves in ^-direction. The interaction 
potentials are g n m{q) for scattering of two incoming elec- 
trons in band n to band m, and h mn {q) for scattering of 
two electrons from band m and n to band n and m with 
momentum transfer q, 

»(«) 



£ 9n 



2e 2 



2e 2 



dydy'K (q\y-y'\)<j> m (y)(f> n (y')(t) m (y')(t) n (y), 

(2a) 

= J dydy'K (q\y-y'\)cj) 2 m (y)ti(y% (2b) 

where K o(z) is the modified Bessel function, e* the effec- 
tive dielectric constant, and 4> n (y) the (real) wave func- 
tions in transversal direction. The vertex corrections to 
the backward scattering are singular when the wave num- 
bers of the incoming and outgoing electrons are close to 
one of the Fermi wave numbers ±fc„, while those to the 
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forward scattering are singular onJw. through the diver- 
gence of the potentials at small qu'H . For the moment, 
we assume short range potentials g nrn and h nm regu- 
lar at q = and show below how to treat the Coulomb 
case. We write the vertices in the form <& s (P\iP2iP3iPa)i 
where p — (pi,£i), Pi, S\ and p 2 , £ 2 are the wave num- 
bers and the Matsubara frequencies of the incoming elec- 
trons, j»3, £3 and pi, £4 are those of the outgoing elec- 
trons, and the suffix s indicates the spins states of the 
incoming electrons, i.e., s = 1 for anti-parallel spins and 
s = 2 for parallel spins. Two kinds of vertices have to be 
distinguished. The vertex 7snm(Pi>P2)P3>P4) ( n 7^ m ) 
is that for p x = -k m , p 2 = k n , p 3 = fc„, p 4 = -k m , 
and x s ™ m (Pi,p 2 ,p 3 ,p 4 ) is that for p x Si -k n , p 2 = k n , 
P3 — k m , Pa — —k m (the case of n = m will be in- 
cluded in x s „ m (p 1 ,p 2 ,P 3 ,P 4 ))- In the parquet diagram 
methodtf , the vertex 7 snm is composed of two irreducible 
parts a snm and S snm : a snm (S smn ) is the sum of the 
all the graphs which can be separated by cutting two 
antiparallel (parallel) electron lines. In the same way 
irreducible vertex parts T snrn and T snrn for Xsnm are in- 
troduced. The corresponding integral equations have the 
form (m ^ n) 

o-i, lm (p 1 ,p 2 ,P3,p 4 ) = ginm 

+ kBT^T{Sxnm(Pi, k,p 3 , k - q)G n {k)G m {k - q) 
k 

X 72nm(fe - 9,P 2 ,fc,P 4 ) +Tl mn (Pl,k,p 3 ,k - q) 

x G m (k)G n {k - q)x2nm(k - g,p 2 ,fc,p 4 ) + (!<-> 2)}, 



(3) 

where q = p 3 p l = p 2 - p 4 , gi nm = g nm (k n + k m ), 
and G n (k) is the one electron Green's function G n (k) = 
1 / (%he v — E n ^k + p), p being the chemical potential. The 
notation (1 «-> 2) means that the terms in which the 
subscript 1 and 2 of the vertices are interchanged are 
to be added. The sum J2k = S e / dk/2n gives rise 
to the logarithmic divergence at zero temperature when 
q = (k n + k m , 0). The parquet graph method takes into 
account the contributions of all graphs giving rise to such 
singularities. The vertices can be regarded as functions 
of two parameters 77 and £, where (v n are the Fermi ve- 
locities) 



f 



V 



log 



r 4 
-0 



(lo 2 + v 2 Aq 2 )(Lo 2 + v 2 n Aq 2 ) 

and a corresponding definition for £ with (Aq,Lo) re- 
placed by (Ap, Ae), where Aq = p 3 — pi — k m — k n , lo = 
£ 3 -£i, Ap = pi+P2 + k m -k n and A£ = £i+£ 2 , e being 
the cut-off frequency of order of Fermi energy devided by 
h. Thus one has cri„ m (p 1 ,p 2 ,p 3 ,p 4 ) = a lnm (^il) etc. 
We introduce A nm := (irh(v n + u m )) _1 and the prod- 
uct F ■ G(C, rf) = Jo F{£, min[t, rj\)G(£, t)dt for £ > 77 and 
F ■ G((, 7]) = $ F(min[i, £] , rj)G(t, rf)dt for r\ > £. In the 
case 77 > £, eq. (3) can be written in the form 

Cl„ m (£, f]) — gi nm — —Kim {Slum ' 72nm(£, r ?) + 

+Ti 

mn ' \2nm 

(e )?? ) + (i^2)}. (5) 

The equations for the remaining vertices are similar ex- 



cept for the spin structures and the terms of lowest or- 
der in the interaction. Since the expressions are quite 
lengthy, we will present the detailed calculations else- 
where. The important point is that the set of integral 
equations can be transformed into differential equations 
whose physical meaning becomes transparent in view of 
the RG calculation below. The result is 

7lnm 

nm L'Tlnm (0 + X 1 nm (0X2nm(£)] i ( 6a ) 
7 2 nm(£) = -Km [lLm(0 + Xlnm(0 + xlnmiO] , ( 6b ) 
XlnmiO = b/lnm(QX2nm(€) + l2nm(0Xlnm(0] 

+ Y, XU iXlnl(0X2lrn(0 + - 2 X l„l(0Xllm{0] . 



(6c) 



X2nm(0 2A nm ^ ' 7s«m K)Xsnm(^) 



s=l,2 



(n 7^ m) 



(6d) 



Xinn(0 = -2^A HX ?„i(0 + 2X) A «XmKC)X2i„(0, 



(6e) 
(6f) 



(5) = -A (C)+^A HX L(0- 
The initial values of the vertex functions are 7inm(0) = 

Xlnm(0) = fflnmi 72nm(0) = 5lnm — ^nmi andx 2nm (0) = 
yinm 92nm- 

We now outline the RG calculation. There, it turns 
out that the case ?j = £, eq. (||), namely Aq = Ap = 0, 
and A£ = u = e\ — £3 = E\ + £ 2 , is sufficient to extract 
the logarithmic divergence of the vertex and to reunp- 
duce eqs. (6). If the energy variables are chosen as intiP, 
£1 = (3/2)o;, — £ 2 = £3 = £ 4 = (l/2)w, the vertex de- 
pends on one single variable lo only (the energy transfer 
in a scattering event) which, however, is sufficient to ex- 
tract the low-temperature dependence from a T = cal- 
culation via Hlu — > kgT. The diagrams contributing to 
7 s „ m and Xsnm to lowest order are the standard bubble, 
triangle, ladder and crossed types. For their evaluation, 
we define Cooper and zero sound functions 

( 4 ) C mn (u>,q):= [ f .„ G m (k,u')G n (q - k,ui - u')u(k), 



(2n) 2 



(7) 



and correspondingly Z nm (oj,q) with (q — k, lu — lo') re- 
placed by {k + <7, 0/ + lo). Here, u(k) is an arbitrary non- 
singular function, and G n (k,Lo) is the free T = 0-Green's 
function. One can easily extract the logarithmic diver- 
gent part for small momentum transfers, C m n{Lo,q — 
k n - km) = -i [u(-k m ) + S nm u(k n )} X nm hx(w/s ). For 
q = k n + k m , one has Z nm (uj,q) = iu(—k n )\ nm lii(ui/s ). 
In both cases, the logarithmic divergence is due to a small 
region in k- space around the points ±fc n ,±fc m , where 
the dispersion E n ^ can be linearized. In the Tomonaga- 
Luttingcr model this linearization is performed from the 
very beginning, here it appears naturally in the cvalua- 
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tion of the diagrams. Careful evaluation yields (n =/= to) 



Xlnm( w ) — 9lnm ~ln \ 2^mn[ginmh'2nm + (1 <-> 2)] 



— A;; [g2nl9llm + .91? 



(8) 



and similar expressions for the other vertices. Here, for 
simplicity we did not include a spin-dependence of the 
interaction parameters or the \ mn as, e.g., in a spin-orbit 
split system or in presence of a magnetic field, although 
the formalism easily incorporates these effects. 

In the renormalization procedure, the cutoff £o pro- 
vides a natural means to scale to the new set of mod- 
elaiiP defined by a smaller value e' < eq. At the same 
time, one requires invariance of the scattering properties 
of the system which means invariance of the vertices, e.g. 
in the case of Xinm 



Xlnm(w, r nm (£o)) — Xinm ( — LO,F nm (e' ) 

\ £ 



(9) 



impurity (backscattering) V nm (k n + k mi T) = V nm (T), 
which becomes temperature dependent and divergent for 
T — > 0, drastically. In the parquet method, upon in- 
serting the full backward scattering vertices, an inte- 
gral equation for V nm {T) is obtained describing the In- 
dependent screening effects on the impurity potential. In 
the equivalent RG calculation, one first starts from the 
perturbative expression for V nm (u>) which again is loga- 
rithmic divergent, so that a scaling procedure must be 
performed, namely 

V nm (u>,r nm (e )) = z(^-)V nm ( (^-)ui,T nm (e' a ) ) . (11) 

Here, scaled and unsealed quantities differ by a multi- 
plicative factor z independent of ui. The bare couplings 
are replaced by the vertices to include all possible inter- 
mediate interaction-induced scattering processes which 
leads to 



where T nm (eo) denotes the coupling matrix with the cor- 
responding 'initial' conditions g S nm{u = £q) = g S nm and 
h snrn (uj = So) = h snm . The functional dependence of the 
vertex on the energy and the coupling constants thus re- 
mains invariant under a scale transformation. For higher 
than second order perturbation theory in the bare cou- 
plings, one should consider the product of the vertex 
with two of its 'arms' (single-particle Greea's-functions) 
as invariant, rather than the vertex itselfclP. However, 
since self-energy corrections £ are divergent only in sec- 
ond order and the bare vertex itself already is first order, 
eq. (^|) is consistent. In fact, the first order Hartree-Fock 
self-energy £ n (fc) merely renormalizes the band ener- 
gies £„fc (counted from the chemical potential), namely 
£,nk = £nk + S„(fc) — /i, with the Fermi momenta k n as 
solutions of £,nk n = 0. 

The dependence of the coupling parameters on e' gives 
the flow to the new set of models with modified cutoff. 
By differentiating eq. (||) with respect to ui and setting 
the new cutoff e' a — to, 

Q 

XlnmM = ( e °/ W )~9Q * lnm ( fi >f( W ))ln= So ' ( 10 ) 

one obtains one of the RG equations eqs. (6), and corre- 
spondingly for 72„ m (w) etc. Frequency dependent singlet- 
and triplet vertices correspond to the respective cou- 
plings in the RG equations, 7i mn (u;) = h lmn (uj) ,72mnM = 

hl mn (u)-h2mn{u),Xlmn(u) = fflmnM, and X2mn(uj) = 

(w) - g2mn(u)- One can easily verify the agreement 
with eqs. (6) upon introducing the dimensionless loga- 
rithmic variable £ = — ln(w/eo)- 

We now turn to the conductivity a of the multi-subband 
wire, namely its temperature dependence, where the power- 
law due to the interaction effects is predicted. The latter 
are expected to be important in relatively clean wires, 
where the impurity potential V nm (q) for scattering from 
band n to m with momentum transfer q can be treated 
as a perturbation. However, the screening by the inter- 
acting electron gas is non-perturbative and changes the 



t— \-nV nm {L0) = . 

ouj to 



(12a) 
(12b) 



CWM := -A nm ^2 pYsmnM + Xsnm(w)] , (12c) 
s=l,2 

for n ^ to, where we defined V nm := V(k n + k m ), and 
"/nm = for the case n = to. For low temperatures or 
correspondingly small w, one can replace the vertices by 
their fixpoints, i.e their £— >oo(w — > 0) values. The solu- 
tion gives the temperature dependence upon replacement 
of uj by ksT/h, 



V nm {T) = V nm -{T/T y 



(13) 



where a nm := a nm (uj = 0) and Tq :— e^/kB- Here, we 
neglected logarithmic corrections of the form log(T/To) 
which experimentally are not very important compared 
to the power-law dependence of the conductivity a to be 
derived now. 

Since the impurity potential is treated perturbatively, 
a can be obtained from the Boltzmann equationEj-* for 
the multi-channel case, 



N 



2e2 

■kTi ^— ' 
n=l 



N 



(14a) 



(14b) 



where LK nm = m* 2 \V nm {T)\ 2 /(k n k m ?f), LJ n 
\Vnm(k n ~k m )\ 2 /(k n k m h 4 ) (to ^ n), and J nn = - ^ m ^„ 
Jnm — J2 m Knm, L being the length of the wire, and 
m* the electron band mass. In the above, we used the 
bare scattering amplitudes for the intersubband forward 
scatterings since their temperature dependences are not 
strong. Below we assume that a nm > for all n and to. 
In fact it is the case for all the examples of our numerical 
calculations discussed below. Then, at very low temper- 
atures the K nm are divergent like K nm oc T~ 2oinm , and 
we may neglect J nm except J nn which contains K nm . 
We obtain ^2 m= iK nm (lm + In) = 1- Note that each 
term in the left hand side of this equation is positive and 
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hence has to remain finite in the limit of zero tempera- 
ture although the K nm are divergent. Let M{n) be the 
value of m for which a nm is the largest for given n, and 
a n = a n M(n) ■ Then, assuming that l n oc T 2a ™ at very 
low T, the term J^m=i -^-nm'n tends to a non-vanishing 
constant in the limit of zero temperature. As for the 
first term J2 m =i K n m l m , we find that K nm = K mn , and 
hence that K nm l m remains finite because K nm diverges 
more slowly than or like T~ 2oim . Thus the assumption 
l n oc T 2q " is reasonable, and the conductivity behaves 
like a oc T 2q , at very low temperatures, where a is the 
smallest of the a n . 

So far we have not taken into account the divergence 
of the bare interaction potential at small q — ► 0, see 
eqs. (2). The elements g n n{q) and h nm (q) are diver- 
gent like h nm {q) = g nn {q) « -(2e 2 /e*) log Wq = U c {q) 
for q —* 0, where W is of the order of the width of 
the wire. This divergence is of the same kind as the 
one which appeared in the self-consistent equations for 
the vertices. Therefore in diagrams containing g n n(q) 
and h nm {q) with small q one has to replace the inter- 
action with the screened Coulomb interaction U sc (q) = 
W(l + ^(9)Eti^(?)) where q = (q,w). Here, 

77,(q) =2 Vl q 2 /[nh(vfq 2 +u% (15) 

is the polarization function in subband I. Typical values 
of lu and q for g n n(q) are those with uj ~ v n q. Since q is 
very small, one can put 

9 N 

On the other hand, the situation for h nm {q) is not so 
simple; we find that the important value of u> is an aver- 
age of v n q and v m q which for h nm (0) leads an expression 
like eq. ( fl6l) with v n q — > (v n + v m )q/2. This is sufficient 
to study the dependence of the exponent a nm on N, 
the number of the occupied subbands, since the problem 
reduces to that of short range interactions and we can 
apply the preceding results. 

In order to obtain the exponents a nm , we have to solve 
the coupled equations eqs. (6). The fixpoints of the latter 
at £ — > oo are obtained by putting the right hand sides to 
zero. The initial values of j snm and Xsnm are real, and 
hence they are real for arbitrary £. Then, from the equa- 
tion for 7 2llm (£) we find that 7i„ m (oo) = Xi»m(°o) = 
X2«m(oo) = 0, (n m) and hence Xinn(oo) = 0. Al- 
though the eqs. (6) are very complex, we can solve them 
in some limit. Considering the expressions eqs. (2), we 
first note that unless Wq <C 1, h nm (q) and g n n(q) are of 
order e 2 /e*. On the other hand, as can be seen from 
cq. (|l6|), 5„„(0) and h nm (0) are of order of the typ- 
ical Fermi velocity times U . Therefore the ratios of 
h nm (0) or g nn (0) to h nm (q) or g nn (q) for q - 1/W are 
of order a* times the typical Fermi wave number, where 
a* = 2e*h/m*e 2 is the effective Bohr radius. The typ- 
ical value of the Fermi wave number is of order l/W, 
and those ratios are typically of order a* /W, as are the 
ratios of X2nn(0) and 72nm(0) to the initial values of the 
other vertex functions. Thus when a* /W ^> 1 (weak 
coupling limit), in the right hand side of eqs. (6) we can 



neglect terms containing neither X2nn(£) nor 72nm(£)- 
^From the remaining nonzero equations we find Xsnrn 
to decrease exponentially if A„„X2rm(0 +A mTO X2mm(£) - 
2A„ m 72rim(£) < 0, where the weak coupling approxima- 
tion is consistent and valid. Then, X2rm(£) and 72nm(£) 
do not depend on £, and with eq. ( |T3| ) it follows a nm = 
KmKmify and a nn = X nn g nn (0). eq. @ yields a nm = 
a nn = 1/2N, and with l n oc T 2q ™ we arrive at our pre- 
diction eq. ([!]) a ~ . We note that the exponent 
1/N no longer depends on the interaction parameters in 
the weak coupling limit. 

In order to verify the validity of the latter, we did 
numerical calculations for the case of two occupied sub- 
bands. The quantum wire was modeled by hard walls 
of width W. We calculated the quantities Xsnm(O) and 
Jsnm(O) and solved eqs. (6) numerically. Detailed results 
will be given elsewhere. We used the values of j2nm(0 
and X2nn(0 f° r £ = 10 instead of the fixpoints 72nm(oo) 
and X2n«(oo)- The Xsnm(£)'s do not depend on £ very 
much up to £ = 10 and we find that the results are at 
large in agreement with a nm — ow — 1/2N. Practi- 
cally, the values of £ are limited by the temperature T or 
by the length L of the wires, £ = 10 corresponds to 10~ 4 
times the typical Fermi temperature. If the latter is of 
order 100 K, the weak coupling approximation is valid 
down to 10 mK. 

Integrals over the jjsave number should be cut off at 
1/1/ in the lower limitils' . This corresponds to a cutoff for 
£ at £l ~ log(l//c„L), k n being the typical Fermi wave 
number. If k n = 10 6 /m ~ 10 5 /m, £l = 10 corresponds 
to L = 100/iTO ~ 10/iTO so that our results are applicable 
to real systems. 
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